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Abstract 

We consider a gravitating system of vanishing cosmological constant consisting of 
an electromagnetic field and a scalar field coupled to the Einstein tensor. A Reissner- 
Nordstrom black hole undergoes a second-order phase transition to a hairy black 
hole of generally anisotropic hair at a certain critical temperature which we compute. 
The no-hair theorem is evaded due to the coupling between the scalar field and the 
Einstein tensor. Within a first order perturbative approach we calculate explicitly the 
properties of a hairy black hole configuration near the critical temperature and show 
that it is energetically favorable over the corresponding Reissner-Nordstrom black 
hole. 
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1 Introduction 

The "no-hair" theorems are powerful tools in studying black hole solutions of Einstein 
gravity coupled with matter. These "no-hair" theorems describe the existence and sta- 
bility of four-dimensional asymptotically flat black holes coupled to an electromagnetic 
field or in vacuum. In the case of a minimally coupled scalar field in asymptotically flat 
spacetime the "no-hair" theorems were proven imposing conditions on the form of the self- 
interaction potential pQ. These theorems were also generalized to non-minimally coupled 
scalar fields [2]. 

For asymptotically flat spacetime, a four-dimensional black hole coupled to a scalar field 
with a zero self-interaction potential is known [3]. However, the scalar field diverges on 
the event horizon and, furthermore, the solution is unstable [4], so there is no violation of 
the "no-hair" theorems. In the case of a positive cosmological constant with a minimally 
coupled scalar field with a self-interaction potential, black hole solutions were found in [5] 
and also a numerical solution was presented in [6], but it was unstable. If the scalar field 
is non-minimally coupled, a solution exists with a quartic self-interaction potential [7] , but 
it was shown to be unstable [HIE]. 

In the case of a negative cosmological constant, stable solutions were found numerically 
for spherical geometries [TO], [TT] and an exact solution in asymptotically AdS space with 
hyperbolic geometry was presented in [12] and generalized later to include charge [13] and 
further generalized to non-conformal solutions [14] . This solution is perturbatively stable 
for negative mass and may develop instabilities for positive mass [15J. The thermodynamics 
of this solution was studied in [T2] where it was shown that there is a second order phase 
transition of the hairy black hole to a pure topological black hole without hair. The 
analytical and numerical calculation of the quasi-normal modes of scalar, electromagnetic 
and tensor perturbations of these black holes confirmed this behaviour [16] . Recently, a new 
exact solution of a charged C-metric conformally coupled to a scalar field was presented 
in [T7J Q2] . A Schwarzschild-AdS black hole in five- dimensions coupled to a scalar field was 
discussed in [TS], while dilatonic black hole solutions with a Gauss-Bonnet term in various 
dimensions were discussed in [20] . 

Recently scalar-tensor theories with nonminimal couplings between derivatives of a 
scalar field and curvature were studied. The most general gravity Lagrangian linear in 
the curvature scalar R, quadratic in the scalar field 0, and containing terms with four 
derivatives was considered in [21J. It was shown that this theory cannot be recast into 
Einsteinian form by a conformal rescaling. It was further shown that without considering 
any effective potential, an effective cosmological constant and then an inflationary phase 
can be generated. 

Subsequently it was found [22] that the equation of motion for the scalar field can be 
reduced to a second-order differential equation when the scalar field is kinetically coupled 
to the Einstein tensor. Then the cosmological evolution of the scalar field coupled to the 
Einstein tensor was considered and it was shown that the universe at early stages has 
a quasi-de Sitter behaviour corresponding to a cosmological constant proportional to the 
inverse of the coupling of the scalar field to the Einstein tensor. These properties of the 
derivative coupling of the scalar field to curvature had triggered the interest of the study 



-3- 



of the cosmological implications of this new type of scalar-tensor theory [23| |2~^| |25| 126] . 
Also local black hole solutions were discussed in [27J . 

The dynamical evolution of a scalar field coupled to the Einstein tensor in the back- 
ground of a Reissner- Nordstrom black hole was studied in [28]. By calculating the quasi- 
normal spectrum of scalar perturbations it was found that for weak coupling of the scalar 
field to the Einstein tensor and for small angular momentum the effective potential outside 
the horizon of the black hole is always positive indicating that the background black hole is 
stable for a weaker coupling. However, for higher angular momentum and as the coupling 
constant gets larger than a critical value, the effective potential develops a negative gap 
near the black hole horizon indicating an instability of the black hole background. 

The previous discussion indicates that the presence of the derivative coupling of a scalar 
field to the Einstein tensor on cosmological or black hole backgrounds generates an effect 
similar to the presence of an effective cosmological constant. In this paper we investigate 
this effect further. We consider a spherically symmetric Reissner-Nordstrom black hole 
and perturb this background by introducing a derivative coupling of a scalar field to the 
Einstein tensor. We show that in this gravitating system there exists a critical temperature 
in which the system undergoes a second-order phase transition to an anisotropic hairy black 
hole configuration and the scalar field is regular on the horizon. This "Einstein hair" is the 
result of evading the no-hair theorem thanks to the presence of the derivative coupling of 
the scalar field to the Einstein tensor. 

The coupled dynamical system of Einstein-Maxwell-Klein-Gordon equations is a highly 
non-linear system of equations for which even a numerical solution appears beyond reach. 
To solve the field equations, we expand the fields around a Reissner-Nordstrom black hole 
solution and pertubatively determine the critical temperature. By solving the first-order 
equations numerically near the critical temperature, we study the behaviour of the hairy 
black hole solution. We calculate the temperature of the new hairy black hole and compare 
it with the corresponding temperature of a Reissner-Nordstrom black hole of the same 
charge. We find that above the critical temperature the Reissner-Nordstrom black hole 
is unstable and by calculating the free energies we show that the new hairy black hole 
configuration is energetically favorable over the corresponding Reissner-Nordstrom black 
hole. 

The paper is organized as follows. In Sec. [2] we set up the field equations and outline our 
solution. In Sec. [3] we find the zeroth order solution and calculate the critical temperature 
near which the Reissner-Nordstrom black hole may become unstable and develop hair. In 
Sec. H] we find the first-order solutions to the system of Einstein-Maxwell-Klein-Gordon 
equations which are hairy black holes near the critical temperature. In Sec. [5] we discuss 
the thermodynamic stability of our first-order hairy solution. In Sec. |6] we discuss the 
validitiy of our perturbative expansion. Finally, in Sec. [7] we conclude. 

2 The field equations 

Consider the Lagrangian density 

C = - ~F, U F^ - {gT + KGnDMDuV)* - m 2 \ V \ 2 , (2.1) 
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where = <9 M — ieA^ and e, m the charge and mass of the scalar field and k the coupling 
of the scalar field to Einstein tensor, of dimension length squared. 

In this paper, we shall concentrate on the case of a massless and chargeless scalar field, 
setting 

m = , e = , (2.2) 

and leave the study of the general case to future work. Consequently, the scalar field ip is 
real. 

With the choice of parameters (12.21) . the field equations resulting from the Lagrangian 
(12. ip are the Einstein equations 

R pv = %nGT^ v , (2.3) 

the Maxwell equations 

VF^ = , F vv = d^A u - d v An , (2.4) 
and the Klein-Gordon equation for the scalar field, 

1 --d, [V^GT + KGT)d v <p] = . (2.5) 



The stress-energy tensor receives three contributions, 

T,u = + T$ + k&, u , (2.6) 

where T^f is the electromagnetic stress-energy tensor, T($ is the standard scalar field 
contribution, 

T$ = d,<pd u¥ > - ^g^d^dpip , (2.7) 

and Q^ u is an extra matter source resulting from the derivative coupling of the scalar field 
to the Einstein tensor, 

= --V^V^ - 1 -{V V ) 2 G^ + iv M V,(Vv?) 2 - i^D(V^) 2 

+^V Q V /3 (V C VV / V) . (2.8) 

The field equations have well-known solutions for ip = 0, the Reissner-Nordstrom black 
holes. We are interested in finding hairy solutions, with (p ^ 0. The no-hair theorem will 
be evaded thanks to the coupling of the scalar field to the Einstein tensor ("Einstein hair"). 

To solve the non-linear field equations, we shall expand around a Reissner-Nordstrom 
black hole. It should be pointed out that not all Reissner-Nordstrom black holes can be 
starting points of such a perturbative expansion. We will show that for a given charge of 
the black hole, there is generally a unique mass for which the Reissner-Nordstrom black 
hole can be a starting point. The corresponding temperature is then expected to be a 
critical temperature at which a phase transition occurs from a Reissner-Nordstrom black 
hole to a hairy one. 
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Introducing the (small) order parameter e, we therefore set 

<p = e (/)+ £ vn...) 

9^ = gf}+e 2 g$ + ... 

A, = Af+e 2 A^ + ... (2.9) 

and solve the field equations perturbatively. 

To find a static solution of the field equations, it is convenient to consider the metric 
ansatz^ 

d s 2 = -e~ a dt 2 + l 2 e a [e-\dr 2 + r 2 d6 2 ) + r 2 sin 2 Odcf) 2 ] , (2.10) 
and electromagnetic potential 

A t = A(r, 6) , A = 6. (2.11) 

The metric functions are expanded as 

I = l + e\ + ... 
a = «o + e 2 ai + . . . 

P = P + e 2 p 1 + ... , (2.12) 

and the electrostatic potential as 

A = A + e 2 A 1 + ... . (2.13) 
The components of the electromagnetic stress-energy tensor are 

rp(EM)t = _rp(EM)<f> = _1_ [r2{drA)2 + 

T (EM)r = _ T (EM)0 = _1_ _ ( ^ )2] 

Ti EM)e = -^d r Ad e A. (2.14) 
I 9 

The Einstein equations reduce to the convenient subset 



R 



\ + R* = SttG (l* + T* - T) , (2.15) 
R\ = SttG/i?- V) , (2.16) 



2 

R r r -R 9 e = 8nG(T;;-TH) , (2.17) 

where 

T = T£ = -{d V f , (2.18) 



1 This is a form of the Lewis-Papapetrou metric written in isotropic coordinates in which all the metric 
functions I, a and /3 depend only on r and 9 |29) . 
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to be solved together with the Maxwell equation (Gauss's Law), 

V^VM = . (2.19) 

and the Klein-Gordon equation (I2.5p . 

The Hawking temperature of the black hole solution can also be expanded as 

T = T + e 2 T 1 + ... . (2.20) 

When e — 0, the temperature is T = To. This is a critical temperature at which a second- 
order phase transition occurs to a hairy black hole, if the latter solution exists and is 
energetically favorable. For a small order parameter, 

£~|T-T | 1/2 , (2.21) 

so our expansion can be viewed as an expansion around the critical temperature. 
Next, we proceed with the perturbative solution of the field equations. 



3 Zeroth order solution 

At zeroth order we obtain the Einstein-Maxwell equations 

< } = SttGT™ , VFJ = 0, (3.1) 
as well as the Klein- Gordon equation 

1 



-g(0) 



y/-g(0)(g^ + nR^ v )8 v p {0) =0, (3.2) 



where we used G^} = R$ for the solution of (13. ip . Thus at this order, the Einstein- 
Maxwell equations decouple (effectively, e = 0), and the solution g$ , has no hair 
(e.g., a Reissner- Nordstrom black hole). 

Note the Klein-Gordon equation for ip^ in the background (g$ , Ay) constrains the 
background, because not all solutions of the Einstein-Maxwell equations lead to regular 
functions (p(°>. The demand of regularity results in relations among the parameters of the 
Reissner-Nordstrom solution. Viewed as a thermodynamic system, the critical temperature 
T is also constrained by demanding regularity of (p(°>. 

The Einstein-Maxwell equations at zeroth order (13. ip have as a solution the Reissner- 
Nordstrom black hole. In Lewis-Papapetrou coordinates, 

u 2 , n io U 2 2ucoth5 
l = 1 - ^ , l e a ° /2 = 1 + ~2 + — , 0o = O, 3.3 

and the electrostatic potential is 



-B 



A? = $ - ^e-^' 2 , Q = -jj± , $ = , (3.4) 

rlo VGsinhB VG 
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where Q is the charge of the black hole, and $ is the value of the potential as r — > oo, fixed 
by the requirement that A^A^ 1 be finite at the horizon. 

To put it in more familiar terms, performing the coordinate transformation 

u 2 

p = l re ao/2 = r + — + 2/i coth B , 
r 

the metric becomes 

ds 2 = -f{p)dt 2 + 4^ + P 2 dQ 2 , f(p) = l-(l+ ^f) ?± + ^f- , (3.5) 



showing that the horizon is at p = p + , 

Pi 



Q^G e 



The electrostatic potential reads A^ = — , confirming that Q is the charge of the hole. Its 
temperature is 

An 47rp + V Pi J 47r/i 
The mass, entropy, and potential are, respectively, 



M^coth*, S^-^- B , (3.7, 



From these, or using the Euclidean action, we deduce the Gibbs free energy 

F = M~TS-Q^ = ^. (3.8) 

Additionally, the scalar field obeys the Klein-Gordon equation f !3.2p . For k — 0, there 
is no regular solution, which is a consequence of the no-hair theorem [TJ [2] . However, for 
k 7^ 0, we obtain regular static solutions for certain values of the parameters \x and B. 
In [28] the Klein-Gordon equation (13. 2p was solved numerically in the background of the 
Reissner- Nordstrom black hole. The corresponding temperature of the Reissner- Nordstrom 
black hole To is the critical temperature, near which the black hole may become unstable 
and develop hair, evading the no-hair theorem. 

To calculate the critical temperature, it is convenient to introduce the coordinate 

z = £ , (3.9) 
r 

so that the boundary is at z = and the horizon at z — 1. 

We look for regular solutions of the Klein-Gordon equation ( 13. 2 j) in the interval [0, 1]. 
Such solutions where discussed in (28] where an instability was found in the quasinormal 
frequencies for multipole number i ^ 0. To capture such an effect we consider a non- 
spherical ansatz 

y? (0) = Z(z)P i (cos9) . (3.10) 
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After some straightforward algebra, we obtain 

z 2 

— (Z (l + kR)Z')' - (1 - kK)£{£ + l)Z = , (3.11) 

H) 

where prime denotes differentiation with respect to z and we defined 

K = RK = = -R^ e = 8sinl f B T? „ / — . (3.12) 

* At 2 [(1 + z 2 ) sinh B + 2z cosh B}* v ; 

The Klein-Gordon equation (13. 2p does not have a solution with a regular scalar field at 
the horizon for £ = 0. This means that scalar hair, if it exists, is anisotropic!! We shall 
concentrate on the case of a dipole, 

£=1. (3.13) 
At the boundary we obtain z 2 Z" — £(£ + 1)Z ~ 0, therefore 

Z ~ z i+1 . 

At the horizon we have ((1 — z)Z')' w 0, therefore 2' ~ or 2' ~ y^. For each value 
of the charge Q, there should be a unique combination of parameters \i and B that gives 
a convergent (and therefore a finite S) at the horizon. For given k, we will determine 
numerically pairs of fi and B from (13.111) and their values will be used to determine T as 
a function of the charge Q using (13.61) . 

To solve the scalar wave equation ( 13. lip numerically, we use the variational method. In 
particular we use an expansion around the horizon 

oo 

Z{z) = z 2 W{z) = z 2 ^W n {l-z) n . (3.14) 

n=0 

It turns out that for a better than 10% precision, one needs to keep at least 10 terms in 
the expansion. 

The coefficients W n depend on k and B. Our aim is to find a unique k which will give 
finite Z(z) and Z(z)' at the horizon z — 1. One may rearrange eq. (13. lip to get 

fjdz(l-2?)[2Z(z)*/z* + Z"(z)] 
/„ dz(l - z 2 )n{z)[2Z{z) 2 /z 2 - Z'\z)\ ' 

which can be solved graphically. For the choice 

k = 10, Q = 9.5, (3.16) 

in units in which 8irG = 1, we obtain /x = 0.15, B = 0.158 (see figured]), which gives the 
critical temperature 

T ps 0.01 . (3.17) 

In figure [21 we depict the scalar function Z(z) for the chosen values of the parameters, nor- 
malized so that Z(l) ~ T (the scalar field and the temperature have the same dimension). 
The normalization is arbitrary, since the wave equation is linear. 
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Figure 2: The function Z(z) for the scalar field versus z for the parameters of figured! 




Figure 3: The critical temperature To vs. the charge Q for various values of k. Black squares 
correspond to k = 15, white squares to k — 10, while there follow 3 and k = 1. 

The function Z(z) is regular in the entire range outside the horizon, as desired. 

For the first order numerical solutions of the next section we used the values of k and Q 
given in ( I3.16p . However, to show how the critical temperature depends on the parameters, 
in figure |3] we show the results of a numerical study of the critical temperature To as a 
function of the black hole charge Q for various values of the Einstein coupling constant k, 
ranging from k — 1 to k — 15. The critical temperature diverges as Q — > (Schwarzschild 
limit) for all values of k. 

4 First order solution 

At first order, the field equations yield a solution near the critical temperature, where 
the scalar field (p backreacts on the metric. 

To extract the first-order equations from the full set of non-linear field equations, notice 
that the stress-energy tensor to this order has contributions from the electromagnetic stress- 
energy tensor (12.141) . and the scalar field (eqs. (12.71) and (12.81) . with ip replaced by ep^ and 
g^ v replaced by the RN solution y$u). For the contribution (12. 8p . after some algebra we 

2 We should point out that spherically symmetric solutions of the Klein-Gordon equation resulting from 
the Lagrangian density (|2.ip may still exist in the case of non- vanishing mass and charge of the scalar field. 
This more general case is under investigation. 
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Figure 4: The function qjio(^) of the angle- independent part of the first-order correction to 
the metric. 



obtain the explicit expression 



ei = e z z = -eg = -e; = e 



2 



2z 6 e 



-3a 



fi% sinh^ B 



z 2 Z' 2 cos 2 6 + Z 2 sin 2 9 



0(e 4 ) . (4.1) 



We give the technical details of the solution of the equations (I2.15p - fl2.19p in Appendix A. 
Having the solutions of the equations ( 12.15p -f j2.19p we will numerically determine the metric 
functions a and and the electric potential A. Notice that the order parameter e and the 
normalization of the zeroth-order scalar field tp^ are not independently defined, since only 
their product (which ought to be small for the perturbative expansion to be valid) enters 
the field equations. For convenience, we set e — 1 in our numerical calculations, making 
the normalization of ip^ small. 

We wand to find the functions ac\a(z), Axq{z), 0:12(2), Ai 2 (z) which appear in the first 
order corrections 



«i = a w (z)P (cos9) + a 12 (z)P 2 (cos 6) , 

Ai = A w (z)P (cos6) + A 12 (z)P 2 (cos6). (4.2) 

We will calculate first the angle-independent first-order corrections. For the correction 
a 10 we work with equation (1A. 12h noticing that at the boundary aio(z) ~ Xz, while at the 
horizon such a solution yields an expression a + 61n(l — z). Fine tuning of A will yield 6 = 0, 
that is, a regular solution. On the technical side we change the variable ctio to £o = — 
with the new boundary conditions Co(0) = A, Co(0) = 0- Figure H] depicts the results for 
a>io for A = —0.000421852, which is found to yield a regular solution. 

The first-order correction to the electric potential A(z) = A + e 2 A w may be directly 
determined if a is known using flA.lip and the result is plotted in figure 
A procedure similar to the one used for the scalar field will be used for the solution of the 
equations (1A.14|) and (1A. 15[) for the angle-dependent first order corrections A\ 2 and ol\ 2 . 



We use the expansions: 

oo oo 

a 12 (^) = z 3 5>iJ(l-zr, A 12 (z) = ^^^(l (4.3) 



n=l n=l 
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Figure 5: The electrostatic potential A(z) at first order. 

where 5^ , A$ depend on a single free parameter which can be chosen to be A^ 2 (it is 
easily seen that the first term vanishes, A^ 2 = 0). This parameter can be determined by a 
variational method. To this end, multiply flA. 14|) by Ayi and integrate. We obtain 



l 

dz 



o 



., , , I 2 2 1 -3z 2 6 \ , 2 Ql . , 



;z 2 L> 2 z 2 L> z 2 Zg ^ 2 y 12 fiD 2 



(4.4) 



where we defined 



D(z) = l + z 2 + 2z coth B . (4.5) 
Similarly, if we multiply flA. 15[) by a 12 , we obtain 







/ dz 


-(«; 2 ) 2 -( 


Jo 





l + ^ 2 6 8 \ 2 16nGQ , cr KQ2 

+ ^ + tt; ■ i mo a i2 H ; — oci2A 12 — IbrrG — — a 2 



/ 2 z 2 (D sinh Bf J " 12 ^ fil ~ 12 ~ 12 XU " ~ ^ 2 







6) 

Both of these equations must be satisfied at the right value of A\ 2 ■ Either one of them 
determines this value. 

In figure [6] we show the results for the variational expressions versus the undetermined 
parameter A 12 and find that A l2 = 0.2. We depict the left-hand side of (14.61) and the sum 
of (14.41) and ( 14.61) keeping four terms in each of the expansions (14. 3D . Notice that the two 
curves are almost indistinguishable and both approach zero at the same value, confirming 
the consistency of our numerical approach. In figure [7] we depict the solution for A12, while 
in figure [8] the solution for an- Therefore, the first order corrections A\i and ol\i are regular 
everywhere. 

Finally, the metric function which is given by flA. 18[) has as the first order angle- 
independent correction the function /3 10 of (IA.23I) which is shown in figure O The function 
Pi in the angle-dependent first-order contribution given by flA. 241) is shown in figure [TOl 
Therefore both corrections are regular in the entire range of z. 



5 Thermo dynamical stability at first order 



In this section we will discuss the thermodynamical stability of our solution. We need 
to know the temperature of the hairy solution. The temperature of the hairy black hole at 



Figure 6: Variational expressions 

dHD and gH) vs. A\1> showing that A\1> = 0.2. 




0.2 0.4 0.6 0.8 1.0 



Figure 7: The function Ai2{z) for the parameters of figure [Hand A 12 = 0.2. 



0.01 - 




Figure 8: The function a 12 (z) for the parameters of figure [Hand A% = 0.2. 
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Figure 10: The function /3 n . 
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first order expressed in terms of the critical temperature is 



Tha . r = Tne£ 2 (-aio(l)+/3io(l)/2) _ 



(5.1) 



Notice that Th air > T , so the RN black hole is unstable at high temperatures (above T ) 
for a fixed charge Q. As the mass M approaches its minimum value at extremality, the RN 
black hole becomes stable. 

Since the correction is quadratic in the scalar field, we deduce for the value of the scalar 
field at the horizon, 

T, 



(p+ = (fi = 7 W— 2^-1. (5.2) 
z=l v J o 

Let us find an RN black hole at this temperature. We shall find one with the same charge 
Q. Call the parameters of this black hole /i' — /i + e 2 fa and B' = B + e 2 B\. Since 
\xj sinh_B = const., we have 

B x =tanh J B^i . 

H 

Another relation between B\ and \i\ is found from setting 



We obtain 



ST 
-«io(l) + 



-aio(l) + 



Mi) 



^o(l) 



2 e 2B - 1 n 
These two relations determine B 1 and fa. For the free energy, we deduce 



SF, 



RN =£2 



RN 



l + e 2B r , m , A>(1) 



Notice that SF RN > as long as e 2B < 3. 

We need to compare it with the free energy of the hairy black hole. The mass of the 
hairy black hole is found from the asymptotic behaviour of a, 



GM 



/i 



«o(0) , 



where prime denotes differentiation with respect to z. Therefore, 

GM hair = GM RN + e 2 ^a' w (0) . 

This is found from the surface (Gibbons-Hawking) term in the action. The entropy is found 
from the area of the horizon, 



S hair = S RN e-^-™ + ^ 
Therefore, the product TS remains unchanged. 



z=l 
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There is an additional contribution from the Einstein-Hilbert action because the Ricci 
scalar does not vanish. We obtain a contribution to the free energy 



8F, 



EH 



- / d 3 xy/^g- 



R 



16nG 2 

This is the dominant change in the free energy and is clearly negative. Explicitly, 

2itu. f 1 dz 



5F, 



EH 



J ^e 2 [{zZ'f + 2Z 2 } < . 



Therefore, the difference in free energies is 

1 



AF — SF hair —6F, 



RN 



1 + e 2B 
-e 2 a' w (0)F RN - 3 _ e2B £ 2 



-aio(l) + 



Ml) 



F RN +6F EH . (5.3) 



Putting numbers in (I5.3P we find AF = —0.08091, showing that the hairy black hole is 
thermodynamically stable. 



6 Discussion of the solution 

To complete the first-order solution and verify the validity of the perturbative expansion, 
we determine the first-order correction to the scalar field and calculate various invariants 
of the metric and show that they are regular at and outside the horizon, showing that no 
singularity arises at this order. 

With our choice of the zeroth-order scalar field as a dipole (eqs. (I3.10p and (13. 13ft ). 
the first-order correction (eq. (12.91) ) contains both a dipole and a i = 3 term. Let 

ip^\z, 9) = (fio(z) cos 9 + tpn{z) cos 3$ . (6.1) 

The field equation obeyed by <pM is obtained by collecting the first-order terms in the 
Klein-Gordon equation (12.51) . The resulting equation is too long to be included here. It is 
straightforward to see that it results into decoupled equations for <p w and <pn. The latter 
involve the functions (p(°>(z), aio(z), au(z), (3w(z) and 0n(z), which have already been 
calculated. Both ipio(z) and fn{z) can be seen to behave as c\qz 2 and c\\z 2 in the limit 
z —7- 0. One may tune the coefficients C\q and Cn to ensure that the corrections vanish 
at z = 1. The results for the two functions are depicted in figure [TTJ It is readily seen, 
upon comparison with figure that the corrections are of the order of the zeroth order 
contribution, so the series in e is expected to have a finite radius of convergence. 

Having demonstrated regularity of the first-order corrections to the scalar field, we now 
turn to addressing the same question regarding the metric. To this end, we need to compute 
gauge-invariant quantities, such as the Ricci scalar. The Ricci scalar R vanishes at zeroth 
order since it corresponds to a Reissner-Nordstrom black hole. We computed R at first 
order both analytically and numerically and found that it is regular everywhere. In figure 
[I2]we plot R versus 9 for representative values of the radius, namely z = 0.3, z = 0.6 and 
z = 0.9. Furthermore, we computed two additional gauge-invariant quantities, R^ V R^ V and 



Figure 11: The first-order corrections of the scalar field, <fio(z) (left panel) and <fn{z) 
(right panel). 

RnvpaR^ vpa at first-order both analytically and numerically and found them to be regular 
everywhere (see figure [13]). We show the values of the product R llv B} lv of the Ricci tensors, 
as well as the product R lllip(7 W 1Vf " T of the Riemann tensors versus 9 for various values of z. 
We find again that the results are finite and roughly of the same order of magnitude for 
the three typical values of z. 

7 Conclusions 

We studied the effect of the presence in the Einstein-Hilbert action of a derivative 
coupling of a scalar field to Einstein tensor, on static black hole solutions. We considered 
the Reissner-Nordstrom black hole solution in isotropic coordinates and in this background 
we introduced a scalar field coupled to Einstein tensor. For small values of the scalar field 
we studied in details how the derivative coupling backreacts on the metric, solving the full 
coupled dynamical system of Einstein-Maxwell-Klein-Gordon equations. 

We found that the Reissner-Nordstrom black hole above a certain critical temperature is 
destabilized to a new hairy black hole configuration. We studied the properties of this new 
hairy black hole solution near the critical temperature and we showed that the scalar field 
is regular on the horizon and at infinity. The no-hair theorem is evaded due to the presence 
of the derivative coupling of the scalar field to the Einstein tensor. This new "Einstein 
hair" solution is in general anisotropic with the scalar field and the metric functions to 
depend also on the angular coordinate. We calculated the mass and the temperature of 
the new hairy black hole solution and by considering the free energies we showed that the 
new hairy black hole configuration is thermodynamically stable. 

It would be interesting to extend the analysis to more general "Einstein hair" by includ- 
ing mass and charge for the scalar hair and explore the existence of spherically symmetric 
hair. In particular, it would be of great interest to see if hair can develop down to zero 
temperature and form a configuration of vanishing entropy. Work in this direction is in 
progress. 
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Figure 13: R^BF" (left panel) and R^ pa R» upa (right panel) vs. 9 for for z = 0.3, 0.6, 0.9 
(bottom to top). 
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A First-order solution of Einstein-Maxwell equations 



In this appendix we give the technical details of solving the field equations fl2.15p - fl2.19p 
at first order in the order parameter e 2 . We start with (12.151) . On the right-hand side, only 
(p(°> contributes and only through Q^ v . Using ( I4.ip . we obtain 



z 2 d 2 z l - zd z l + d 2 l + 2 cot Odgl = 8vrG^^e a (-0* - Gt) + 0(e 4 ) = C(e 4 ) , (A.l) 



which has clS cL solution 

I = 1 - z 2 + 0{e A ) , (A.2) 

therefore I is unchanged at first order, l\ = (recall eqs. (13. 3p and (13. 9p for an RN black 
hole). 

Next we solve the Maxwell equation which takes the explicit form 



z 2 d 2 z A + z 2 d z (a + In l)d z A + d 2 A + [cot 6 + d e a + d e In I] d e A = . 
At first order, 



(A.3) 



z'diA + 



-— + o z a 



z 2 d z A + d 2 e A + cot 6d e A = 0{e 4 ) . 



From (I2.16p . at first order we have 



2z 3 

z 2 d 2 a - —d z a + d 2 e a + cot 6d a = 8irGe a 

i: 



I 2 

2. J Ct 



[zd z Af + 2n 2 K-6 t 



We expand the first-order corrections in Legendre polynomials, 



(A.4) 



on = aio(z)Po(cos9) + cti2(z)P2(cos6 
A 1 = A 10 (z)P (cos6) + A 12 (z)P 2 (cost 



fi 2 ^e a Ql = e 2 [Qo(z)P (cos6) + Q 2 (z)P 2 (cos6)] + 0(s 4 



where 



Qo 



K 
~6 



zZ') 2 + Z 2 



Q2 = j [{zz'f - z 2 ] , 



and 71 is given in (I3.12p . 



(A.5) 
(A.6) 
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Using 
we obtain 



P:'(cos6) + cot 6 P'A cos ( 



(A + e 2 A w ) 



2z 



2z , \i 
— + (a + sa 10 ) 



+ l)P £ (cos( 



(A + e 2 A 10 )' = 0(e 4 



A'' + 



2z 



Otr 



A' - —A, 



-a' A' 



(A.7) 
(A.8) 



(a + e 2 a l0 )" - — (a + £ 2 aio)' - 87rGe ao+£ Ql0 (A + e 2 A 10 ) ,z 



£ 2 Qo 



12 



l " : 
From ( 1A.7I) we deduce 



^a' 12 - — a 12 - 87iGe ao A' 2 a 12 - 167iGe a ° A' A' 2 



z 2 

(A.9) 

167tGk% . (A.10) 



A (z) + e 2 A w (z) 



Q f 1 dz' 



l(z') 



-ao(z')-e 2 a 10 (z') 



(A.ll) 



Eq. ( 1A.9|) implies that the angle-independent part of first-order correction of a satisfies 



2z 



a 



10 



l-z 



a in — 



2 "io 



[(l + z 2 )sinh J B + 2zcosh J B] 



7«io = 16/tG 



(A.12) 



To solve it, first we look at the boundary conditions. At the boundary, a" ~ 0, therefore 

a io ~ \z . (A. 13) 

At the horizon, (1 — z)a'{ — a' w ~ 0, therefore aio ~const., or aio ~ ln(l — z). The solution 
that asymptotes to Xz as z — > yields a mixture ccio ~ a + 61n(l — z) at the horizon. There 
is a unique value of A for which b = and the solution is regular. 

The remaining two equations form a coupled system to be solved for the angle- dependent 
first-order contributions A 12 and a,\i. Explicitly, 



A" 

^12 



'0 



Z 



1 + z 2 + 2z coth B lr 



+ 



A' l2 - —A X2 

Z A 



"12 



Q 



If, 



H {1 + z 2 + 2z coth B) 2 " 12 ' 

(A.14) 

16 ^^ 12 + 16^. 

z 2 

(A.15) 



[(1 + 2 2 )sinh J B + 22cosh J B] 2 



As 2 — > 0, we have 



12 



Az 3 + 



a X2 = &z + 



(A.16) 



H l 



The constants A and a are fixed by demanding a 12 = A 12 = at the horizon (so there is 
no angular dependence of the temperature or an electric field along the horizon). 
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Finally, we obtain (5 from the third Einstein equation f)2.17p . At zeroth order, we have 



= — (za'n 



3^ + 2 



In 



^ + 87tGe ao (zA' ) 2 
In 



(A.17) 



which is easily seen to be satisfied. 

To solve the equation at first order, we set 



and obtain 
where 



Sn(z) 



(1 + z 2 )zP[ = S (z) , 

Az 2 [(l + z 2 ) coshS + 2,2 sinh B] 
(1 + z 2 ) sinh B + 2z cosh B 

%z 2 {\ - z 2 ) 

[( 1 + 2 2 ) sinh B + 2^ cosh 



(A.18) 
(A.19) 



a 



10 



-a 



12 



8ttG(1 - z 2 



Akz 4 



[i 2 [(1 + z 2 ) sinh B + 2z cosh B\ 



to be solved for Pin(z), and 



(l+z i )zP' n -2(l-z 2 )p ll = S 1 (z) 



(A.20) 



(A.21) 



where 



S 1 (z) 



12z 2 [(l + z 2 ) cosh B + 2zsmhB] , 



as, + 



12z 2 (l - z 2 



(l + ^ 2 )sinh J B + 2zcosh J B 12 ' [(1 + z 2 ) sinh 5 + 2zcosh J B] 2 



"12 



+ 8ttG(1 - z 2 )[(zZ'Y + Z 2 ] + 8vrG 



4/^ 4 (l - 



/i 2 [(l + z 2 ) sinh S + 2z cosh B] 



(A.22) 



to be solved for pu(z)- 

They are both first-order equations and are easily integrated. We obtain 



Pw{z) 



dy 



o 1/(1+ y 2 ) 



So(y) 



(A.23) 



where we used the boundary condition /?io(0) = 0. 
For p u , we need Pu = at both ends. We obtain 



P 



ii 



(l + z 



2\2 



(A.24) 



where the limit of integration was chosen so that Pu = at the horizon (z = 1). Notice 
that Pu ~ 0(z 2 ) as z — >• (r — )■ oo), so the other boundary condition is also satisfied. 
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